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The quantum derivatives of e~ A ,A~ 1 and log A, which play a basic role in quantum 
statistical physics, are derived and their convergence is proven for an unbounded positive 
operator A in a Hilbert space. Using the quantum analysis based on these quantum deriva- 
tives, a basic equation for the entropy operator in nonequilibrium systems is derived, and 
Zubarev's theory is extended to infinite order with respect to a perturbation. Using the 
first-order term of this general perturbational expansion of the entropy operator, Kubo's lin- 
ear response is rederived and expressed in terms of the inner derivation 5-h for the relevant 
Hamiltonian TL. Some remarks on the conductivity <j(u>) are given. 
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^ ■ §1. Introduction 

o ' 

Recently, the present author 1 ) ~ 3 ) proposed a new scheme of quantum calculus, 
the so-called quantum analysis. In this scheme, the derivative of an operator-valued 
function with respect to the relevant operator itself is expressed only in terms of the 
original operator and its inner derivation (i.e., a hyperoperator or superoperator), 
Qh' and an operator expansion formula is derived. 

r~| ■ In the present paper, the quantum derivatives of e~ A ,A _1 and log A, which are 

basic operator functions in physics, are derived, and their convergence is proven in 
§ 2 for the unbounded positive operator A in a Hilbert space. (See also Appendices 
A and B.) Nonlinear responses in equilibrium are expressed in terms of quantum 
.^J | derivatives in § 3. A basic equation for nonequilibrium systems is derived in § 4 

^ ■ using quantum analysis. This derivation has the merit that it is valid even for an 

unbounded entropy operator. On the other hand, Zubarev's derivation is based on 
the power series expansion of the density matrix with respect to the entropy operator, 
and consequently it is restricted to a bounded entropy operator. Zubarev's theory 4 ) 
is extended to infinite order in § 5. This gives a renormalized perturbation theory 
with respect to an external field. Kubo's formula of linear response 5 )' 6 - 1 is then 
rederived and expressed in terms of an inner derivation in § 6. Some remarks on 
the conductivity cr(co) are given in § 7. The entropy operator rj{t) in a dissipative 
system [namely — log p(t) for the density matrix p(t)] is expressed in a compact form 
using the inner derivation in Appendix C. This expression is convenient for studying 
quantum effects, because it is expressed only in terms of commutators. 
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§2. Quantum derivatives of e A , A 1 and log^4 and the convergence of 

the differential df(A) 

In a previous paper, 1 ) quantum analysis was formulated in a Banach space, 
namely for bounded operators. The term 'quantum analysis' refers to noncommuta- 
tive differential calculus in terms of inner derivations, namely commutators. Formal 
expressions and several formulas of quantum analysis are derived in Ref. 1). In prac- 
tical applications, for example, to quantum and statistical physics, we often have to 
treat unbounded operators in a Hilbert space. As is well known, it is difficult to prove 
generally the convergence of such formal expressions for unbounded operators. 7 )' 8 ) 
Fortunately, the density matrix p in statistical mechanics is a contraction operator 
when the relevant Hamiltonian fi is unbounded (even for a finite system) but posi- 
tive definite (or bounded below). Furthermore, a perturbation may often be assumed 
to be bounded in statistical physics. (For example, a Zeeman energy is expressed 
by a bounded operator in a finite system, while the kinetic energy of an itinerant 
electron system is unbounded.) Hereafter we discuss the quantum calculus for these 
situations. Thus we study here the convergence of the Gateaux differential 

= lim HA + MA)- HA) ^ 
h^o h 

where A is unbounded but f(A) is bounded, and dA is an arbitrary bounded operator 
independent of A. To consider this situation is one of the key points for studying 
the convergence of Eq. (2-1). The quantum derivative df(A)/dA is defined 1 ) by 

df(A) = ^dA. (2-2a) 

Here df(A)/dA is a hyperoperator which is a function of both A and the inner 
derivation 5 a defined by Eq. (2-7b). This property is crutial in quantum analysis. ^ 
In fact, we have the formula 1 ) -3 ) 

m. = 6 -m (2 . 2b) 

dA 5 a 

in a Banach space. Higher-order quantum derivatives will be discussed in § 3. 

i) Quntum derivatives of e~ A and A^ 1 

Here, we attempt to prove the convergence of Eq. (2-1) for two typical operator 
functions, f(A) = e~ A and f(A) = A -1 , where A is a positive (but unbounded) 
operator. Clearly we have 

|| e~ A ||< 1 and || A' 1 ||< oo, (2-3) 

under the condition that A > a > for a constant a. 
First note that 7) ' 9) 

± e -(A+xB) = _ f 1 e ~(l-s)(A +X B) Be - S (A +X B) d ^ (2 . 4&) 

dx Jo 
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Integrating Eq. (2-4a) we obtain 

e -(A+xB) = e -A _ f X dt f l dse - {1 - s)iA+ tB) Be - S (A + tB)_ (2 . 4b) 
JO Jo 

Then we can prove the convergence 

lim || {e' A+hB) -e- A )/h- f 1 e- ( - 1 -^ A (-B)e- sA ds \\= (2-5) 

when A is a positive (but unbounded) operator and B = dA is bounded, as is shown 
in detail in Appendix A. Thus we arrive at the differential 

d(e~ A ) = - f 1 e-^~ s)A (dA)e~ sA ds = -e~ A A(A)dA, (2-6a) 
Jo 

or the quantum derivative 

dp~ A r 1 

— = -e- A + J o e-^ A 5 cxp{ _ sA) ds = -e~ A A(A). (2-6b) 

This is well defined for a positive operator A. Here, the hyperoperator A(A) is 
defined by 1 ) 



A{A) = f 1 e t5A dt = - - - , (2-7a) 
Jo 



e 5A - 1 

with the inner derivation 8a defined by 



5 A Q = [A, Q]=AQ- QA. (2-7b) 

The ratio of the hyperoperators (e 5 - 4 — 1) and 5a is well defined, although d^ 1 does 
not necessarily exist. The formula (2-6b) with Eq. (2-7a) will be used frequently 
later. 

Concerning the convergence of the power series expansion of e~ A A(xA), we have 
the theorem. 

Theorem 1 : The power series expansion ofe~ A A(xA)dA with respect to x converges 
in the uniform norm topology for A > and for | x |< a' 1 , where a is defined by 
the upper limit 

a= lim || (A^SATdA ||£ . (2-8) 

n— >oo 

The proof is easily given using the Stirling formula n\ ~ ri n e~ n for large n and 
the following inequality. 

Inequality : When A > 0, we have 

|| e~ A 5 r XB || < n n e~ n \\ {A~ l 5 A ) n B \\ (2-9a) 
for any positive integer n. 
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The proof of the above inequality is easily given using the inequalities 

II e~ A S%B || < || e~ A A n || • || (i" 1 ^)^ || (2-9b) 

and 

|| e~ A A n || < e~ n n n . (2-9c) 
It should be noted here that 5 A and A commute. 

Corollary : If B x l k is defined for any positive integer k and there exists the max- 
imum M = maxfc || A~ 1 B 1 / k A ||, then the power series expansion of e~ A A(xA)B 
with respect to x converges in the uniform norm topology for A > and for 
| x |< 1/(M + 1). 

Proof: First note that 

{A- l 5 A ) n B = ]T(-l) fc ( n j A~ k BA k (2-10a) 
for any positive integer n. Then we have 



||< ^ (™ j || A~ k BA k 



fc=0 

n 



< 2 (?) H A- l B l l k A \\ k < (M + 1)™ (2-10b) 
k=o \ k J 

under the conditions of the above corollary. Thus we obtain 

lim || (A- l 5 A ) n B \\k < M+l. (2-10c) 

n— >oo 

Note that lim^oo || A~ 1 B 1 / k A || = 1. Then the maximum number M may exist 
when the deformation of B 1 ^ by the transformation of an unbounded operator A is 
finite for all values of k. 

Similarly we study the differential of the resolvent operator A^ 1 when A > a > 0. 
We easily obtain 



lim || ( A \ „ - I") //» - 

h-o " \A + kB A) 1 A y 'A 



< lim | h | • || -r || 2 • || „ 1 , - || • || B || 2 = 0. (2-lla) 



That is, we have 



d(\) = = {-A' 2 + A-H A -,)dA = - A{A l _ 8A) dA. (2-llb) 

This gives 

d .1, 1 , 
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This is also bounded when A > a > 0. Here we have used the relation 1 ) S^-i = 
A- 1 -(A -5a)- 1 . 

ii) Quantum derivative o/log^4 

The above arguments can be extended to more general case in which f(A) is 
also unbounded but df(A) is bounded for bounded B = dA. A typical case is given 
by f{A) = log A. The operator log(^4 + hB) is formally expressed by the following 
integral 



= / 7 7 )dt + h / -B 

Jo \t + l t + A Jo t + A t + 



—dt 
A 



~h 2 I -^—B-^—B } , dt. (2-12) 

■ t + A t + A t + A + hB v ; 



Then we have 



-dt 
A 



ll |logM + ftB) _ loM1/ , ! _^_L_ B _L 

roo 1 1 

<\h\-\\B\\ 2 II— ~r II 2 -II 7—r-z\\<tt- ( 2 - 13 ) 

~ Jo 11 t + A 11 "t + A + hB" x ' 

Consequently we arrive at 10 ) 

Clearly this is bounded when A is positive (i.e., A > a > 0) and dA is bounded. 
This is formally written as 

^ZT = 3 " f TTA W-* = -tf Ml - (2-15) 
The second expression of Eq. (2-15) gives the convergence of df(A)/dA. 

iii) Convergence of df(A) for an unbounded operator A and for the bounded differ- 
ential dA. 

In general, the derivative df(A)/dA is formally given by the following formula. 1 ) 
Formula 1 : When f{x) is an analytic function of x, we have 
df(A) _ 5 f(A) _ f(A) - f(A - S A ) _ f 1 , (1 



f fW(A-t6 A )dt. (2-16) 
Jo 



dA 5a 5a 
Here f^ n \x) denotes the rath derivative of f(x). This is formally expanded as 

^ = f {l \A) - y^\A)5 A + ■■■ + 7^J {n+l) (A)5l + ■ ■. (2-17) 



6 



Then, we have the following theorem. 

Theorem 2 : Let A be unbounded, and let {f( n \A)} for n = 0, 1, 2, • • • and dA be 
bounded. Then, the formal expansion (2-17) operating on dA converges to Eq. (2-16) 
in the uniform norm topology if 

f {n+l XA) i 

ai= lim II i- \f-S%dA ||»< 1. (2-18a) 

n + 1! A 11 y ' 



n— >oo 



A proof of this theorem is easily obtained. Theorem 1 is a typical example 
of the above general theorem. This theorem can also be extended to higher-order 
derivatives (see § 3 and Appendix B). 

In the more general situation in which the operators {f^ n+1 \A)5^dA} are un- 
bounded, the convergence proof of Eq. (2-17) can be studied using the strong norm 
convergence. Then, the condition (2-18a) is replaced by 

a[ = lim II L \>s n A dA^j \\n<\ (2- 18b) 

1 n-+oo 11 (n + 1)! 
for ip € V with some appropriate domain V in Hilbert space. 

§3. Higher-order quantum derivatives and nonlinear responses in 

equilibrium 

The higher-order quantum derivative d n f{A)/d n A is formally expressed ^ by 
the multiple integral 

= n\ J o dh ^ dt 2 ■ ■ ■ J Q dt n fW(A - t 1 S 1 t n 5 n ), (3-1) 

where f(x) denotes the n-th derivative of f(x) and the inner derivation 5j is defined 
by 

Sj-.dA-dA dA = dA ■ dA {5 A dA) dA. (3-2) 

Then we have the following operator Taylor expansion formula : 1 )' 2 ) 

f(A + xB) = ±^^m.. B n (3 . 3 ) 
n=0 

with the notation B n = B B . 

It is sometimes important to study nonlinear responses in condensed matter 
physics, as in spin glasses (in which only nonlinear susceptibilities diverge 11 )' 12 ) at 
the transition point). 

As is well known, an equilibrium system is described by the canonical density 
matrix 

p = e-W- 11 ® (3-4) 

for the Hamiltonian H of the system in the presence of an external field H conjugate 
to a physical quantity Q. When Q does not commute with H, nonlinear responses are 
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described in terms of the canonical correlations of Q, namely by a multiple integral 
of the time correlation function of Q using the Feynman formula. They are now 
expressed as 

00 (-H) n rPp-PK 

n=0 

in quantum analysis. Thus the n-th order nonlinear response is expressed by the 
n-th order quantum derivative of p. The above expression (3T) of higher derivatives 
of f(H) = e~@ n in terms of the inner derivation 5% is convenient for evaluating 
the required nonlinear responses explicitly, for example, using the high-temperature 
expansion method. The above static perturbational expansion with respect to the 
external field H can be extended to that of the nonequilibrium density matrix p{t) 
given by a solution of the von Neumann equation (4-1). 

§4. Basic equations in nonequilibrium systems 

As is well known, the density matrix p(t) in a nonequilibrium system satisfies 
the von Neumann equation 

ihj t p(t) = m),p(t)]=6 H{t) p(t) (4-1) 

for the time-dependent Hamiltonian 7i(t) of the relevant system. 
Now we attempt to find a solution of the exponential form 

p(t)=e- r >®. (4-2) 

Concerning the "entropy operator" 77(f), we have the following formula which was 
pointed out by Zubarev. 4 ** 

Formula 2 : The entropy operator 77(f) defined in Eq. (4-2) satisfies the equation 

ih^ = [H(t), V (t)}. (4-3) 
This is a simple example of the following general formula. 3 ) 

Formula 3 : Any operator- valued function f(p(t)) of the density matrix p(t) satisfies 
the equation 

ihj t f(p(t)) = [H(t),f(p(t))}- (4-4) 
It is instructive to give here a compact proof due to quantum analysis: 

dt Jyiy " dp(t) dt dp(t) 

= ~^dpW 6 ^ m = - 5 /W*»W(*) = Mt),f(p(t))]. (4-5) 

Here we have used Eq.(2-16). The above equation (4-3) is our starting point for 
deriving the renormalized expansion scheme (5-10). 
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§5. General perturbation theory on the entropy operator in 

nonequilibrium systems 

We formulate here a general perturbation expansion of the entropy operator for 
the Hamiltonian Ti.it) taking the form 

H(t) = H-AF(t), (5-1) 

with a time-dependent external force F(t) (as in Kubo's linear response theory 5 ) ,6 )). 
Here A denotes an operator conjugate to the external force F(t). Now, we define 
the correction term r/(i) in 

Ti(t)=$ + PH + rf(t) (5-2) 
for r)(t) = — log p(t), where (3 = l/k^T and <P is a normalization constant such that 

e^ = Tr e -/™-V(-°°). (5.3) 

We then expand the correction term rf(i) as 

00 

r,\t) = J2rin(t), (5-4) 

71=1 

so that r] n (t) is of nth order in F(t). This is a new type of renormalized perturbation 
theory for nonequilibrium systems, because even the first-order term r]i(t) gives 
partially infinite-order terms in the density matrix p{t). It is easily shown from 
Formula 2 that rf(t) satisfies the inhomogeneous equation 

j t v'(t) = ±[H(tU(t)}-(3F(t)A (5-5) 
with the initial condition rj'(—oc) = 0, which corresponds to the condition 

p(-oo) = p eq = e-^/Tr e~ m . (5-6) 
Here we have also used the relation 

A = ^-[A,n] = ^s A n. (5-7) 

in in 

Equation (5-5) is the basic formula derived here. A new aspect of this equation is 
that it has the temperature-dependent source term — f3F(t)A. Since H(t) contains an 
external force F(t), Eq. (5-5) is nonlinear with respect to this force. The linearized 
equation is given by 

j t vi(t) = ^[n, m (t)]-pF(t)A. (5-8) 

The solution of Eq. (5-8) with the initial condition r/i (—00) = is obtained as 
771 (t) = -13 y* F(s)exp Q^{t - s)5n^j Ads = -(3 J° e £S F(t + s)A{s)ds. (5-9) 
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The adiabatic factor e £S has been inserted to insure convergence. 

The above first-order approximation {pi(t) = exp[— <P — 0H — iji(t)]} gives 
Zubarev's statistical operator 4 -* when the Hamiltonian 7i(t) is given by Ti(t) = 
TL — AF(t). This first-order approximation, namely Zubarev's theory, is justified 
if the second-order term 772 (i) is much smaller than r]i(t). 

For higher-order correction terms of r/(t), we have the following. 

Formula 4 : The higher-order entropy operators {r] n {t)} are given by 

m(t) = ~^ [° e £S dsF(t + s) f F(t + s')[A(s'),A(s)]ds' 
in J -00 Jo 

f S f* 1 ftn-2 

Ut) = - W ^ 1 ]_J s dsF { t + s)j od t 1 l dt 2 ...j o dt n ^ 

xF (t + t 1 ).-. F(t + in-i)^( tl )^(t 2 ) • • • dA^^Ais) (5-10) 

with the hyperoperator S^t) an d with A = (ih)~ 1 5A'H- 

These formulas can be derived from Eq. (5-5). They will be useful in studying 
nonlinear responses, because they are renormalized perturbational expansions in 
contrast to the ordinary perturbational expansion 5 ) of the density matrix itself. In 
fact, even the above p\{t) includes terms up to infinite order in F{t). Thus, our 
formulation (5-10) is a new useful result, compared with the ordinary expansion 
scheme of p(t) itself. The quantum analysis of dissipative systems 13 )' 14 ) will be 
presented in Appendix C, using ordered exponentials and free Lie elements. 15 )' 16 ) 

§6. Linear response in terms of the inner derivation 

In this section we discuss linear response as an application of the general per- 
turbation theory presented in the preceding section, and we express it in terms of 
the inner derivation 5% for the relevant Hamiltonian TL. 

The density matrix p(t) for the Hamiltonian H(t) in Eq (5-1) is given by 

up to first-order of in external force F(t). Here, r)i(t) is given by Eq. (5-9). The 
quantum derivative de _/3W / 'd{0H) is expressed by 

im=- e ~ mAimh (6 ' 2) 

as is seen from Eq. (2-6b) . Thus, the first-order term Ap(t) is given by 

Ap{t) = -e^e-^OSWM*). (6-3) 
The average of the relevant current operator J = A is expressed as 
(J) t = TrA P (t)J = -((A(0H) Vl (t))J) 
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= p[° e £S F(t + s)(A(0H)J(s))J)ds 

J — oo 
roc 

= e- £S F{t-s){{A((3H)J(-s))J)ds 
Jo 

= Re{a{u)Fe iwt ) (6-4) 



under the assumption that Tr Jexp(— (3H) = and F{t) = Fcos(wt). Here < • • • > 
denotes the average with respect to the equilibrium density matrix, and the general 
conductivity cr(uj) is expressed as 



rOO 

a{u) = (3 / e~ £S ~ iuJS {(A((3TC)J)J(s))ds; (6-5) 
Jo 



namely 



a( „ ) = m A( m j) —± Jm J } 

= £]T<(/A(/?H)J)[(— SnTJ}) (6-6) 

n=0 

for the Planck constant h, using the hyperoperator A(A) defined by Eq. (2-7a). 
It is also interesting to note that we have 

/ e xn Je~ xn d\= / e xs ™d\J = (3A((3H)J (6-7) 
Jo Jo 

for the current operator J in our notation. This may be thoght of as the 'dressed 
current operator', due to quantum fluctuation. Thus, Kubo's canonical correlation 
(J : J(t)} is expressed as 

(J : J(t)) = - f\e xn Je- xn J(t))d\ = {(A(j3H)J)J(t)). (6-8) 
P Jo 

Then, the Kubo formula for the frequency-dependent conductivity <r(u) is expressed 
in the form 

roo i 

a(u>) =0 (J: J(t))e-^dt = 0{(A{0H)J)-. ttj-^-J). (6-9) 

In particular, we obtain 

a(u) ~ ^-{(A(pH)J)J) (6-10) 

for large to. Some remarks on applications of Eqs. (6-6) and (6-9) will be given in 
the succeeding section. 

The present derivation of the Kubo formula may be more transparent and the 
algebraic structure that <t(oj) is expressed only in terms of the commutators of H 
and J (namely free Lie elements) is convenient in practical calculations, as will be 
shown elsewhere. 
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§7. Some remarks on the conductivity a(uj) 

It is instructive to give some remarks on applications of the formulas (6-6) and 
(6-9) for the conductivity a{u). 

When the current J is a constant of motion, <r(0) is infinite 17 )' 18 ) as seen from 
(6-5). We consider the following more general situation that the current operator J 
contains some (not necessarily all) constants of motion {Hj}, that is, 

J = J2 a J H j + J', (7-1) 

j 

where J' is defined by the remaining part of J orthogonal to all the {Hj}; namely 
J' is off-diagonal with respect to energy 18 ) (i.e., (m | J 1 \ n) = for E m = E n with 
the energy eigenvalues {E n } of the Hamiltonian TL even in a degenerate case). Here, 
the coefficients {aj} in (7-1), namely the ergodicity constants, are given 18 ) by 

«j = {JHj)l{H% (7-2) 
using the orthogonality condition 

(HjH k ) = (H])S jk . (7-3) 
Thus, the zero frequency (or static isolated) conductivity defined by 

roo 

o-(O) =13 (J: J{t))dt 
Jo 

is seen to diverge as 

roo 

(7(0) = pJ2 °? J (Hj)dt + (finite) -» oo, (7-5) 

when at least one of the ergodicity constants {aj} is non- vanishing. This remark 
is useful in practical applications 6 )' 19 ) of the Kubo formula to some exactly soluble 
systems 17 )' 18 ) with an infinite number of constants of motion. 

§8. Concluding remarks 

The quantum analysis introduced in previous papers 1 ) ~ 3 ) has been extended to 
the case of an unbounded operator A in a Hilbert space by restricting our consider- 
ation to the three typical operator functions e~ A , 1/A and log^4 under the situation 
that the differential dA is bounded. The proof is rather easy but it is instructive for 
studying more difficult cases for unbounded operators. 

Our new expressions of response functions in terms of the inner derivation 5n 
(or the dressed current operator A{(3TL)J) are convenient for analytic and numerical 
calculations of these response functions. This result should be compared with the 
abstract operator representation of a KMS-state by Naudts, Verbeure and Weder 20 ) 
in the more complicated situation of infinite systems. 



(7-4) 
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The renormalized perturbation scheme of the density matrix pit) is one of the 
new results in the present paper. This is in sharp contrast to Kubo's well-known 
systematic expansion formula 5 ) of p(t) itself, rather than log pit). 

It is also interesting to note that the quantum analysis is useful in expressing 
an exponential product of a dissipative density matrix in terms of a single expo- 
nential (namely the generalized BCH formula) composed only of commutators, as is 
exemplified in Appendix C. 

Transport coefficients are also expressed in terms of commutators of the relevant 
current operators. 
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Appendix A 

Convergence of Eq. (2-5) 



The convergence of Eq. (2-5) is shown as follows: 

|| {e -(A+hB) _ e -A )/h _ r 1 e - {1 - s )A(_ B)e -sA ds || 

JO 

- [''(It j'' ds L-(ls)(A+tB) B{e - S (A+tB) _ Q -sA } 

h Jo Jo L 



<ll - / '// / </ 

+ { e -(±-s)(A+tB) _ e -(l- s )Ax Be -sA 

|| h || • || B || 2 r rl 

< max 

2 \t\<\h\ 



o 



X 



ds || e -(±s)(A+tB) || 
d\ || e -( s ~ X KA+tB) || . || e -X(A+tB) || 

(A-l) 



+ ['ds f 1 "dX || e -(ls-WA+tB) || . || e -X(A+tB) e -sA 

Jo Jo 

when A is a positive (but unbounded) operator and B = dA is bounded. Therefore, 
we arrive at Eq. (2-5). 

Appendix B 

Expansion Formulas and Convergence of Higher- Order Derivatives 



The rath derivative of f(A) is given by Eq (3-1), namely by the following integral: ^ 
= n! j o d tl ... j o dt n f^(A - (B-1) 
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Here, Sj is a hyperoperator defined by Eq. (3-2), namely by 

Sj : {dA) n = (dA) j ' 1 (5 A dA)(dA) n - j . (B-2) 

This is also formally expanded as follows. 
Formula A : 



d n f(A) 



£ = t ^^f^ m \A) jT 1 dtl • • • /J"" 1 ^ n (E tA-r- (b-3) 



cL4 

m=0 



For example, the first derivative df(A)/dA is given by Eq. (2-17), and 



m=0 v 



2 



d 3 /(^) ^ 3 !(-ir (m+3) 5T +2 



^ 3 ^ (m + 3)! J v ^^(^ + 5 3 ) 
_ (6i + 5 2 r+ 2 ( Sl + 6 2 + 6 3 ) m + 2 

62S3 (62 + $3)63 h ' 1 J 

The expansion of 

^/(^) = = (dA) n (B-5) 



converges in the uniform norm topology when 

a„= lim || ^ ±-L dh--- dt n (Y t i 5 j ) m (dA) n ||-< 1. (B-6) 

Appendix C 

Quantum Analysis of Dissipative Density Matrices 

It is instructive to discuss first the non-dissipative unitary case, 
(i) Unitary case. Here we discuss the von Neumann equation 



ihj t p(t) = [H(t),p(t)}=8 n{t) p(t) (C-l) 

for the time-dependent Hamiltonian H(t) of the relevant system, as in (4-1). A 
formal solution of Eq. (C-l) is given by 

p(t) = exp + ^ J 5 H{s) ds^j p(0) 

- exp + (s Si n(s)ds ) " (0)exp - ("S jf w(s)<is ) ' (C ' 2) 

Here, we have used the following ordered exponentials : 13 )' 14 ) 

ft ft ft ftn-l 

exp + / A(s)ds = l+ / A(s)ds-\ h / dh--- dt n A(t\) ■ ■ ■ A(t n ) H 

JO JO Jo Jo 

(C-3) 
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and 

ft rt rt rtn-1 

exp_ / A(s)ds = 1+1 A(s)ds H h / dh • • • / dt n A(t n ) ■ ■ ■ A(ti) H . 

Jo Jo Jo Jo 

(C-4) 

Clearly the ordered exponentials 

exp + {kfo n{s)ds ) and exp - {-hIo n{s)ds ) (c,5) 

are both unitary and consequently are bounded when {Ti(s)} are self-adjoint. Thus, 
the arguments in § 2 can also be applied to these ordered exponentials, namely 
operator functional. 21 ) It is then shown that the operator functional derivation 21 ) 
for the variation 6H(ti) 



dF[H{t)] tl EE d 



exp + 



1 I* 

^ / H(8)dt 

in Jo 



SF[H(t)] 



tl 5H(h) 

= exp + ^ jf W(s)ds) <5H(ti)exp + ^ jf * W(s)ds^ (C-6) 



is bounded when the elements of {T~t(t)} are self-adjoint and the elements of {5H.(ti)} 

j_ 

ih 



are bounded. Similarly the operator functional derivation of exp_[— ^ J t Tl(s)ds] is 



bounded under the same conditions. 

(ii) Dissipative case. We discuss here the unnormalized density operator p(t) of a 
dissipative system described by the master equation 



dp{t) 



dt ih 



hn,p{t)] + Ap{t) + p{t)A^. (c-7) 



Here, A and A* denote some bounded operators expressing a dissipative effect. The 
Hamiltonian 7i may be unbounded. The normalized density matrix pit) is given by 
p{t) = N(t)p(t) with N^y 1 = Tr p(t). A formal solution of Eq. (C-7) is given as 
follows. First we put 

p(t)=exp^w)/(t)e Xp (-^w). (C-8) 
Then, Eq. (C-7) can be rewritten as 



am 

dt 

where 



= A,f(t) + f(t)4, (C-9) 



• 4 ' = exp (-S W )' lexp (s W ) (C - 10) 
and A\ = (A t )^. Next we put 

/(*) = exp + (- J* A\ds^ g(t)exp_ Qf ' A\ds^j (C-11) 
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with g(0) = /(O) = p(O). Then Eq. (C-9) can again be rewritten as 

dg(t) 



dt 



C(t)g{t), 



where 



C(t) = exp_ {^J A\dsj (A t + A\)exp + J A\dt 
A solution of Eq. (C-12) is given by 

g(t) =exp + [J* /te) g(0). 



Thus we arrive at 



where 



p(t) = exp + C(s,t)ds^j p(t,0), 



£(s, t) = exp (J^ij exp + (- j A\ds ) C{s) 
A\ds ] exp | — — 7i 



and 



x exp_ [ j A' s dsJ exp ^ 



p(t, 0) = exp (^ft) ex P+ (- Aids ) p(0) 
x exp_ ^ yljtis^ exp ^-^-H j . 



(C-12) 
(C-13) 

(C-14) 
(C-15) 



(C-16) 



When both p(0) and A are bounded, p(t) is also bounded. 
Now we put 

p(t, 0) =e-"(*' ). 

Then, we have 

p(t) = exp + ^ £(s,t)ds^ exp(-7/(t,0)). 
Our purpose here is to find the logarithm of Eq. (C-19). For this, we put 



exp + ^ C(s,t)dsj exp(— rj(t, 0)) 



a *(x) 



(C-17) 

(C-18) 
(C-19) 

(C-20) 



Clearly we have $(0) = — rj(t, 0). By differentiating Eq. (C-20) with respect to x, we 
obtain 

e*(*)2\(-<P(x))^^ = £(x,t)e^ x \ 



dx 



(C-21) 



This is transformed into the equation 



d${x) 
dx 



= A- l (-<P{x))e- d *^C{x,t) 



A- l {${x))C{x,t) 

-£{x,t)- & 



<P(x) nl __ ^ loge d *M 



e°*(^) - 1 



e°*(^) - 1 



C(x,i), 



(C-22) 
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using the identity 5& = loge 5 *. Then we can apply Eq. (C-20) to Eq. (C-21). Thus 
we finally arrive at the following formula. 

Formula B : The entropy operator fj(t) of the system described by Eq.(C.7) is 
expressed in the form 

fj(t) = -0(t) = - log p(t) 

log[exp + Cg S C ( s ,t)ds)exp(-5 ri (t fi ) ] 
exp+Uo s C(s,t)ds)exp(~S vm ) - 1 ' 
The final expression (C-23) is much more convenient than 



V(t, 0) - f d ^^f^p^^ C^ t). (C-23) 
Jo exp , J n 5c(s,t)ds)exp(-5 ri(m ) - 1 



fj(t) = - log 



exp + C(s, t)ds^J exp(— rj(t, 0)) 



(C-24) 



because Eq. (C-23) is expressed in terms of the commutators of {C(s, t)} and rj(t, 0), 
namely free Lie elements 15 -* . 

The present formulation can be easily extended to the following more general 
dissipative system : 



dp(t) 



±-[H(t),p(t)} + Mt)p(t) + p(t)A\t). (C-25) 



dt %h 

The expression (C-23) is convenient for studying quantum effects 16 ) in non-equilibrium 
systems. 
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